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Abstract — The linear theory of chiral surface waves in 
cylindrical containers with an elliptical cross-section is 
studied in detail. General solutions for phase and the free 
amplitude are given in terms of Mathieu functions. Our 
numerical results show the dependence of the numberwaves 
and natural frecpiencies on the chiral media with 
eccentricity £. The well known case of a circular media for 
8 = 0 is retrieved and remarkable crossings of the mode 
frequencies for certain values of 8 are found. The 
frequency shift is evaluated numerically. 
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I. INTRODUCTION 

The problems of membrane oscillation have been the 
subject of many papers [1-4]. The theoretical and numerical 
investigations presented in these papers concern 
rectangular, circular and arbitrary shaped media. In these 
cases, non-homogeneous problems can be solved by using 
Green’s function method. The problem of oscillation of a 
non-uniform or non-regular shaped system is solved by 
using approximate (numerical) methods [5]. The method of 
fundamental solutions can be used as an example of such 
approximate methods [4]. For chiral media of regular 
shapes (rectangle, circle), an exact solution to 
electromagnetic problems can be derived. 

The application of Green’s function method to the 
oscillation problem of a chiral media, which occupies a 
finite region in the plane, requires the knowledge of Green’s 
function for the Helmholtz equation. The method of 
fundamental solutions applied to electromagnetic problems 
of arbitrary shaped also requires knowledge of the 
fundamental solution of the Helmholtz equation in the 
plane. Green’s functions of the Helmholtz equation in 
regular regions are well known. These functions for 
problems in rectangular and circular regions with various 
boundary conditions are given in [6-8]. In this paper, the 
derivation of Mathieu’ s function for the chiral Helmholtz 
equation in an elliptical region is presented. In order to 
separate the variables, elliptical coordinates are introduced. 


The solution is obtained in the form of a series of eigen 
functions of the associated boundary problem. 


II. CHIRAL ELECTROMAGNETIC 
FORMULATION 

The Maxwell equations for the macroscopic free 
electromagnetic fields, (without charge and current) are well 
known. We often write Maxwell's equations in terms of 
electric and magnetic fields, E and B. 

VxE = - — B , VB = 0 

3t 


VxH 


— D+J, V-D = p 

at 


(la) 


(lb) 

These equations, however, are not complete. Six more 
equations, the constitutive relations, have to be added 
relating the electric field E , the magnetic induction B , the 
displacement field D and the magnetic field H to each 
other. These constitutive relations are completely 
independent of the Maxwell equations. The Maxwell 
equations involve only the fields and their sources. The 
constitutive relations, however, are concerned with the 
equations of motion of the constituents of the medium in an 
electromagnetic field [9-14], 

We often write Maxwell's equations in terms of electric and 
magnetic fields, E and B. 

defined by with the non locality definitions of Bom- 
Fedorov, [9-11]: 

B =|i c [H + T c (VxH)] 

( 2 ) 

D =e c [E + T c (VxE)] 


(3) 


where T c is the chiral factor. The wave propagation is 

(1 - (0 2 |4 c e c T c2 ) v x (V x E)-2T c [i c e c co 2 V x E-co 2 ([i c e c )E=0 
(4) 
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The particular case with 1 — OTq.g.T 2 = 0 corresponds to 
e:b with E=i H so we have 

2T L VxE-E=0 

(5) 

which satisfy the chiral homogeneous Helmholtz equation 

(v’e)+(H__) 2 e = o 


(5’) 


For each component of E or H , is valid the equation 

(v 2 <j>)+(-!-) 2 4>=o 

2T C 

(6) 

This result corresponds to self-dual solutions to the 
Maxwell equations. These solution are known as 
“instantons,” have gained recognition among experts in 
gauge field theory and mathematical physics. The condition 
that a field configuration is self-dual is not invariant under 
the parity transformation r — » -r because of the opposite 
parity properties of the electric and magnetic field; the 
mirror-image configuration is anti-self-dual. As will 
become clear, the physically relevant configurations are 
represented by a sum of self-dual and anti-self-dual 
solutions, which is invariant under the parity 

transformation. In this paper we discuss this type of 
configuration in elliptical geometry in two dimensions. 

III. PROBLEM FORMULATION IN 
ELLIPTICAL COORDINATES 

To solve the problem, we introduce elliptical coordinates in. 
the interior of the domain D (x, 0) , which are coupled 
with Cartesian coordinates by the following relationships 
(Fig. 1): 



Fig.l: Introduce elliptical coordinates (%, 0) from 
Cartesian coordinates. 
x = dcosh%cos0 


y = dsinh%sin0 


(7) 


where X > 0 , 0 < 0 < 271 □ . The equation of ellipse (7) in 
the elliptical coordinates is: % — % 0 

i b 

, where % 0 = artgh — . The Laplace operator in elliptical 
a 

coordinates has the form 

1 


V 2 U = - 


3 2 u 3 2 ur 


d (cosh" x~ cos" 0) 30" 

( 8 ) 

Differential equation (1) and boundary condition (2) in the 
elliptical coordinates 
are as follows 


1 


3 2 3 2 


1 


d 2 (cosh 2 %-cos 2 0) ^3% 2 30 2 ^ ^2T c ^ 

(9) 

U(x o ,0,T c ) = 0, for O<0<2;t 

( 10 ) 

The solution of boundary problem (9), (10) can be found in 
the series form 

u(x. e. T) = XL I ■*>„ <x. ejr„ (T. ) 

n=l 

(if) 

where functions <I> |rm (X,0) satisfy the homogeneous 
Helmholtz equation 

+4^™ (X> 0) +(-^) 2 d 2 (cosh 2 %-cos 2 0)d> n „ (x, ( 
dX 30" 2T^ m 

( 12 ) 

In our case we are looking for solutions of Eq. (12) in the 
interior of the domain D with boundary condition 

* m (Xo.e. TC ) = 0. for O<0<27rand 0 < X ^ Xo 

(13) 

In order to derive eigenfunctions < T > mn (X, 0) , the method 
of separation of variables will be used. We assume that 

0c> 0) = nrjx, q m „ )v(e, q mn ) 

(14) 

where N is a normalization constant. After the substitution 
of (11) into equation (9) and separation of the variables, two 
equations are obtained: 
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^ !L -(a-2q nm cosh2x)R m =0 

d%~ 

(15) 

,2 

- ( a - 2q inn cos 29)¥ m = 0 

d0 

(16) 

where a is the separation constant and 

(1 / 2T^ n ) 2 = 2q mn / d 2 . From physical considerations 
the function \J/(0) must be periodic, i.e. 

\j/(0) = \)/(0 + 27t) .Taking into account boundary 
condition (13) and equation (14), we found that q mn are the 
roots of the equation 
R m(Mm„) = 0 

(17) 

Moreover, we assume that functions ¥ m ( 0 >q in „) are 

periodic with period □ □ or 2 J . This property holds for 
particular values of separation constant a , which depends on 
the values of q . 

Equations (15) and (16) are well known as the radial 
Mathieu equation and the angular Mathieu equation 
respectively. The pairs of the independent solutions of these 
equations are radial and angular Mathieu functions [2]: 


R m (% 5 qmn) 


jCe m (X,q mn ) 

lse m+1 (X,q nm ) 


(18a) 

where Ce m and Se m+1 Sem are the modified Mathieu 
functions of first kind or radial solutions. 


Vm(0.qnm) 


fce m (e,q mn ) 

lse m+1 (B,q mn ) 


(18b) 

m = 0,1 , 2 

Using (11) and introducing functions 

me 2m(0>q) =Ce m( 0 >q)’ me 2m + l(0>q)= se m + l(0>q)’ 


m = 0,l,2.. (19) 


Me 2m (X,q) = Ce m (X,q), 

Me 2m + l(%>q)= Se m +1 (x.q) 

( 20 ) 

function ^nm(X.0) can be written in the form 


^mn (X. 9) = Me m (X> q„m ) me m (0. ' 0 nm ) 

( 21 ) 

The angular Mathieu functions create the set of the 
orthogonal system in interval [0,2 □], i.e. the following 
orthogonality condition holds: 

1 2n 

- j me m (0,q)rne n (0,q)d0 = 5 mn 

n o 

( 22 ) 

where 8^ Dis the Kronecker delta. This leads to the 
statement that eigenfunctions < S ) mn (%, 0) given by (11) 
satisfy the condition 

JJ (X, (X, 0 )d%d0 = TtS^ 

D 

□ □□□ 

Following the spatial symmetry of Eq. (12), the complete 
set of solutions can be chosen as even or odd 
eigenfunctions. 


IV. EIGENVALUES IN ELLIPTICAL 
COORDINATES 

In elliptic coordinates the condition 
d R m (X>qmn) / 8n = 0 at (x, y) in 2- D reduces to 

d R m (X,q mn )/d%U o =0 

(24) 

with % = arctan h(b / a) , to work with number wave, let 
(1 / 2 r P )d = dk m =£k m , £ being the eccentricity and 
k m = ak m = a / ( 2T^ ) , then from Eq. (24) we obtain the 

— — even 

set of even and odd eigenvalues for k m .n, where k m ,n and 

— odd 

km,n are solutions of 

— even — odd 

Ce m (Xo’£km, n / 2) = 0 , Se m (x 0 ,£k n ,n/2) = 0 

(25) 

with n = 1, 2, . . . . Recall that the eccentricity e of an ellipse 
is related to its semi-axes by the expression 

£ = Vl-(b/a) 2 so that one can write 

X 0 = arc tanh \J\ — £ 2 , and therefore, 

the characteristic values k m , n depends solely on the 
eccentricity £ . Thus, if the ratio value between semi-axes 
b/a is modified, the pattern of the wave amplitude will be 
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different carrying the information of the surface symmetry 

and thereby of elliptic geometry. 

The eigenvalue problem (25) states that for a given m the 
symmetry (odd or even) remains valid for any eccentricity 

value. Thus, eigenvalues k m ,n for even and odd solutions 
with the same m but different elliptical oscillation number n 
cannot cross (anticrossing) as a function of the geometry 
factor 8 , i.e. the symmetry is preserved. Moreover, 

eigensolutions belonging to different Hilbert subspace can 
cross. These important facts can be clearly 
seen in Figs. 2 and 3. Fig. 2 shows the variation of the first 

10 eigenvalues k m ,n as a function of 8. Notice that in the 

limit 8 — > 0 the degeneracy of the modes, that corresponds 
to basins with circular cross sections, is clearly displayed, 

i.e., km,n = a / (2T^ n ) tends to (a / (2T^ n )) with(e = 0) 

An important conclusion arises from Fig. 2 where it can be 

seen that k m ,n increases as the eccentricity increases, i.e. 
narrowed ellipses or chiral surface waves more confined, 

lead to larger values of the eigenvalues k m ,n . Moreover and 
according to the general properties of the Mathieu 

— odd 

functions, k m ,n are more sensitive to the confinement than 

— even — odd — even 

the even eigenvalues k m,n , and the inequality k m,n > k m,n 
is valid for any eccentricity value 



£ 


Fig. 2: Eigenvalues k m ,n as a function of the eccentricity e. 
Dashed (solid) lines represent eigenvalues for the odd 
( even ) modes. 
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Fig. 3: Mode frequencies as a function of the eccentricity 8 . 

A comparison of the elliptic case with the classical circular 
basins assuming the same cross-section area shows a 
strongly changing behavior of the natural frequency G) m 
as a function of the geometric factor 8 . 

Fig. 3 is devoted to the evaluation of C0 m for the first 10 
modes as a function of the geometrical parameter 8 . The 
behavior of (0 m versus 8 is very diverse. As a general 

trend, the even modes (m >0, n = 1) show decreasing 
frequencies as 8 increases. Also, it is observed modes 

reaching to a minimum frequency value for certain £ min , 
increasing (Q m for 8 > 8 min , while there are modes 
exhibiting monotone increasing behavior of (Q m as the 

eccentricity increases (in Fig. 3 the odd modes (1, 1), (1, 2), 
(2, 2) and the even mode (0, 1)). Notice that for the odd 
mode (3, 1), G) mn is almost independent of 8 . A more 

interesting feature is the accidental degeneracy obtained at 
certain values of 8 . From the figure the crossing points of 
two different mode frequencies for large values of the 
eccentricity are clearly observed. 

V. CONCLUSIONS 

The derivation of the Mathieu function of the chiral wave 
equation in an elliptical region with the Dirichlet boundary 
condition has been presented. In order to solve the problem, 
elliptical coordinates were introduced. The function has the 
form of a double series of Mathieu functions, which are 
eigenfunctions of the Helmholtz operator in the considered 
elliptical region. Although, the solution is obtained for the 
Dirichlet condition at the boundary ellipse. 

The solution to the non-homogenous problem of the chiral 
oscillation can be presented in an exact form. The linear 
theory of chiral surface waves in cylindrical containers with 
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an elliptical cross-section was studied. General solutions for 
phase and the free amplitude are given in terms of Mathieu 
functions. Our numerical results show the dependence of 
the natural frequencies on the chiral media and the 
eccentricity 8 of the container cross-section. The well 
known case of a circular media for 8 = 0 is retrieved and 
remarkable crossings of the mode frequencies for certain 
values of 8 are found. The frequency shift is evaluated 
numerically. 
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